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AN  ANALYSIS  OF  AXIAL  ACOUSTIC  WAVES 
IN  A  COLD-FLOW  ROCKET 


F.  E.  C.  Culick 

and 

G.  L.  Dehority 
Research  Department 


ABSTRACT.  The  presence  of  waves  in  a  rocket  chamber  depeqcls 
on  the  relative  importance  of  energy  supplied  by  interaction 
with  the  combustion  processes  and  energy  lost  through  dissi¬ 
pative  processes.  As  an  aid  to  assessing  the  latter,  the 
cold-flow  rocket  appears  to  be  a  promising  tool  for  both 
research  and  development.  Combustion  is  absent  and  acoustic 
waves  are  excited  mechanically  in  a  chamber  through  which  a 
steady  flow  is  maintained.  Hie  analysis  presented  here  is 
intended  to  serve  as  a  basis  for  interpreting  such  measure¬ 
ments.  Only  the  case  of  axial,  or  longitudinal  waves  is 
treated,  and  results  show  mainly  how  the  losses  at  the  exhausi 
end  of  the  chamber  may  be  inferred  from  measurements  of  the 
transient  response  and  the  frequency  response  of  a  chamber. 
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NOMENCLATURE 


A  Admittance  at  boundary,  defined  by  Eq.  9 

a  Sonic  velocity,  also  "effective"  plane  where  mean  chamber 
conditions  are  reached  (Fig.  l) 

b  "Effective"  plane  of  nozzle  entrance  (Fig.  l) 

0^  Specific  heat  at  constant  pressure 

D  Dimensionless  expression  defined  by  Eq.  35 

e  Unit  axial  vector 
z 

f  Dimensional  frequency;  also  "non- homogeneous"  part  of  pressure 
boundary  equation  (Eq.  8) 

h  "Non- homogeneous "  part  of  pressure  wave  equation  (Eq.  7) 

ML 

k  Non-dimensional  wave  number  for  steady  oscillations  -~ 
kg.  Non-dimensional  wave  number  for  purely  axial  mode  oscillations 
kjg  Non-dimensional  wave  number  for  classical  mode  oscillations 
L  Chamber  length 
l  Axial  mode  number  (k^/n) 

M  Mach  number 

Mq/  Magnitude  of  Mach  number  perturbation  at  head  end  of  chamber 
N  Dimensionless  expression  defined  by  Eq. 

£  Normal  unit  vector 
p  Pressure 

Q  Resonance  quality  factor 


R  Radius  of  chamber 


r  Dimensionless  radius  (r/k) 

S  Cross-sectional  area  of  chamber 
c 

S  Nozzle  entrance  area 
P 

s  Dimensionless  entropy  (7|  -p'/p) 

t  Time;  non-dimensional  following  Eq.  2  L (real  time)/(L/a)l 
u  Velocity 

X  Square  of  the  absolute  value  of  complex  non-dimensional  nozzle 
admittance  +02) 

x  Dimensional  axial  distance  from  head  end  of  chamber 
Z  Dimensionless  head-end  impedance  (T|/Mq#) 

z  Dimensionless  length  measured  from  head  end  of  chamber  (x/l) 

a  Attenuation  constant,  also  dimensionless  group  in  Eq.  43  and  real 
part  of  4  (Eq.  45) 

S  Dimensionless  group  in  Eq.  43,  also  imaginary  part  of  4  (Eq.  45) 
Y  Ratio  of  specific  heats 

5  Shift  from  closed  tube  resonant  frequency  due  to  the  presence  of 
a  nozzle  (fi^-/Tr) 

e  Shift  from  resonant  frequency  at  half -power  point  Op-"*)] 

Q  Non-dimensional  space  coordinate  in  a  moving  coordinate  system 

(z  -  M  t) 

'  c 

T|  Non-dimensional  pressure  (pVvp) 

9  Imaginary  part  of  non-dimensional  nozzle  admittance 
A  Imaginary  part  of  k 
\  Wavelength 

U  Real  part  of  non-dimensional  nozzle  admittance 
p  Density 

t  Non-dimensional  time  (used  in  variable  transformation) 
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i^tajewasaaifia^^  iSataS^aev  e^«!i9MK*tt 


4  Phase  angle  between  velocity  and  pressure  perturbations,  also 
tan-^h(-An),  see  Eq.  29 

Classical  closed  chamber  pressure  mode  for  purely  axial  oscillations 
♦K  Classical  closed  chambf pressure  mode 
0  Real  part  of  k 

ft  Non-dimensional  frequency  defined  as  0/(1  -  M)2 

C 

(u  Dimensional  angular  frequency 

SUBSCRIPTS 

c  Chamber  property 
l  Pure  axial  mode 
max  Maximum  value 

N  Classical  closed  chamber  mode 
n  Nozzle  condition 
o  Head- end  condition 
p  Peak  value 

r  Real  part  of  complex  quantity 
+  Wave  traveling  in  positive  direction 

-  Wave  traveling  in  negative  direction 
i  Half-power  point 

SUPERSCRIPTS 

(l)  Imaginary  part  of  complex  expression 
(r)  Real  part  of  cooqplex  expression 
*  Perturbation  quantity 

-  Mean  value 
Vector  quantity 


vi 
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I.  INTRODUCTION 


An  understanding  of  the  behavior  of  waves  in  a  rocket  chamber  rests 
essentially  on  knowledge  of  the  various  sources  and  sinks  of  energy. 

The  primary  source  of  energy  is  of  course  the  propellant  itself  and  the 
driving  mechanism  is  associated  with  a  coupling  between  the  wave  motions 
and  the  combustion  process.  It  appears  possible  to  isolate  this  aspect 
of  the  problem,  so  far  as  solid  propellant  rockets  are  concerned,  not 
only  analytically,  but  also  experimentally  (Ref.  1  and  2).  With  suf¬ 
ficient  development,  the  T-bumer  and  L*-burner  may  provide  adequate 
means  for  studying  the  coupling  (interpreted  often  as  an  "admittance 
function")  without  firing  a  complete  rocket. 

But  the  source  of  energy  is  only  half  of  the  problem.  It  is  always 
present,  and  in  order  that  a  rocket  burn  stably,  there  must  exist  ways 
for  wave  energy  to  be  dissipated  at  a  rate  greater  than  that  at  which  the 
combustion  process  provides  energy.  In  some  respects,  it  is  theoretically 
more  difficult  to  assess  the  losses  than  the  gains.  Indeed,  it  is  possi¬ 
ble  to  compute  much  of  the  energy  losses  only  for  relatively  simple  geome¬ 
tries,  and  even  in  those  cases,  approximations  are  necessary.  Hence  one 
would  like  also  to  be  able  to  measure  the  losses  separately. 

Since  the  losses  associated  with  a  chamber  are  strongly  influenced 
by  the  geometry,  which  in  turn  affects  the  mean  flow,  it  is  in  fact 
necessary  to  work  with  a  chamber,  perhaps  scaled  down.  The  problem  is 
simplified  by  using  an  external  source  of  cold  gas,  rather  than  burning, 
for  the  flow  in  the  rocket.  Acoustic  oscillations  may  be  excited,  also 
by  external  means,  and  the  attenuation  of  the  waves  inferred  in  a  manner 
described  below.  Perhaps  the  s implest  example  of  this  kind  of  device  is 
illustrated  in  Fig.  1  (see  Ref.  3  and  4).  The  chamber  is  fitted  with 'a 
nozzle  which  normally  is  choked.  Air  flows  through  the  porous  plate  at 
the  head  end,  and  oscillations  are  excited  by  a  rotary  valve  upstream  of 
the  porous  plate.  An  alternative  way  is  to  drive  the  oscillations  with 
a  small  speaker. 

This  configuration  evidently  simulates  an  end-burning  rocket,  which 
seems  quite  restrictive.  However,  the  measurements  will  provide  infor¬ 
mation  principally  about  the  damping  of  waves  due  to  the  action  of  the 
exhaust  end  of  the  chamber,  i.e.,  the  effect  of  the  nozzle  and  also  non- 
uniformities  in  the  flow  just  upstream  of  the  nozzle  entrance.  It  is  not 
possible  to  separate  these  in  the  experiments.  Clearly  the  results  are 
applicable  to  other  kinds  of  rockets  having  similar  exhaust  configura¬ 
tions  .  The  more  serious  limitation  is  that  in  the  measurements  reported 


1 


In  Ref.  3  and  4,  only  axial  modes  have  been  treated.  In  this  report  as 
well,  only  axial  modes  will  be  discussed,  but  much  of  the  treatment  can 
be  extended  to  radial  and  tangential  modes. 


FIG.  1.  Sketch  of  a  Cold-Flow  Rocket. 


The  following  remarks  are  concerned  with  an  analysis  of  the  axial 
modes,  with  a  purely  axial  mean  flow.  It  Is  not  the  intent  here  to  com¬ 
pute  the  actual  damping  to  be  expected,  for  in  fact  this  can  be  done  only 
for  simple  cases.  Rather,  it  is  more  important  to  determine  how  the 
formal  analysis  can  be  vised  as  an  aid  to  interpreting  the  experimented, 
data,  thereby  yielding  useful  numerical  results  which  cannot  be  obtained 
by  calculation  alone . 

First,  a  relatively  general  treatment  of  the  problem  will  be  covered, 
based  on  earlier  work  (Ref.  5);  this  leeuls  to  formulas  for  the  frequency 
and  attenuation  coefficient  for  freely  decaying  waves.  Then  traveling 
and  standing  waves  in  a  uniform  mean  flow  will  be  treated.  The  analysis 
Is  specifically  for  the  case  when  waves  are  driven  at  one  end  of  the 
chamber.  Formulas  for  the  frequency  response  and  the  standing  wave  ratio 
are  found.  Perhaps  the  most  important  result  is  the  relation  between  the 
attenuation  coefficient  and  the  properties  of  the  frequency  response. 

The  conditions  under  which  the  familiar  definition  of  "Q"  is  valid  are 
found. 

The  three  calculations  mentioned  in  the  previous  paragraph  correspond 
to  the  three  kinds  of  measurements  most  easily  made:  direct  measurement 
of  the  attenuation  of  waves,  the  frequency  response,  and  amplitude  of 


of  pressure  as  a  function  of  position  in  the  chamber  (i.e.,  standing  wave 
ratio).  Thus  the  results  are  immediately  applicable  to  experimental  work. 


II.  GOVERNING  EQUATIONS  FOR  ACOUSTIC  OSCILLATIONS 
IN  COLD-FLOW  ROCKETS 


A  rather  general  analysis  of  acoustic  oscillations  in  rocket  chambers 
exists  (Ref.  5)  and  it  will  be  specialized  in  this  section  to  the  simpler 
case  of  a  cold-flow  rocket  with  axial  mean  flow.  Only  axial  modes  are 
discussed,  but  it  is  necessary,  for  a  reason  apparent  shortly,  to  retain 
the  equations  in  three-dimensional  form.  Extension  to  other  inodes  may 
be  carried  out  by  following  the  analysis  of  Ref.  5-  It  is  convenient 
for  the  analysis  of  the  acoustics  of  rocket  chambers  to  use  an  equation 
for  the  pressure  rather  than  temperature,  so  that  the  three  equations 
for  inviscid  flow  of  a  perfect  gas  are 

■^  +  V  •  (pu)  =  0 

TET  +  U  *  +  2E.  s  0 

Ot  p 


+  YpV  •  u  +  u  •  Vp  =  0 


The  last  may  be  found  by  combining  the  usual  energy  equation  with  the 
continuity  equation.  In  the  usual  way,  all  quantities  are  written  as 
sums  of  mean  values,  independent  of  time  but  in  general  functions  of 
position,  and  fluctuations  which  vary  with  both  time  and  position  in  the 
chamber.  The  complete  perturbation  equations  may  be  put  in  the  form 


where  a  is  the  mean  speed  of  sound,  T)  =  p  Vvp,  and  s '  =  T]  -  p'/p  is  the 
fluctuation  of  entropy.  For  applications  to  the  low  Mach  number  flow  in 
a  cold  rocket,  it  is  not  necessary  to  include  variations  in  the  mean 
pressure  and  temperature .  The  equations  have  the  simpler  form 
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+  ft  •  w'  +  5'  •  ^ 

.  -4  -4 

■^  +  M  •  ^T)  +  Y'i'P  •  ft  +  9  •  M/  «  0  ( 4 ) 

All  variables  are  nov  dimensionless^  t  Is  the  real  time  divided  by  L/a, 

L  being  the  length  of  the  chamber,  ft  «  u/a,  M'  «  uVa,  and  all  lengths 

are  normalized  with  respect  to  L.  Hereafter,  V  •  M  =  0  will  be  assumed. 
This  is  a  consequence  of  the  low  speeds  Involved,  but  it  cannot  be 
assumed  if  there  is  combustion  In  the  chaaber. 


III.  APPROXIMATE  SOLUTION  TOR  THE  FREQUENCY  AND 
ATTENUATION  OF  STANDING  WAVES 


The  results  obtained  In  this  section  will  apply  to  waves  which  are 
not  driven  by  external  means .  A  wave  equation  for  T)  Is  found  in  the 
usual  way  by  taking  the  time  derivative  of  Sq.  4  and  substituting  'Eq.  5 

for  dyi'/dt 

^T)  -  .  M  •  V  -  V  •  [m  •  VM'  +  M'  •  \7ft] 

dt 

The  boundary  conditions  on  "H  are  set  In  accord  with  Eq.  J 
VT|  =  -  -  [m  •  VM'  +  M'  .  w] 


For  ease  in  writing,  the  last  two  equations  will  be  written.  In  the  case 
of  steady  waves,  as 


where 


V2!)  +  k2Tl  »  h 

n  •  VT|  x.  -f  (on  z  =  0,  z  =  l) 

h  =  ikft  •  VT)  -  V  •  •  VM'  +  *'  • 

f  *  lkATI  +  (ft  •  VM'  +  i?'  •  •  n 


(5) 

(6) 

(7) 

(8) 
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The  unknown  boundary  values  of  the  Mach  number  (or  velocity)  fluctuations 
have  been  replaced  in  Eq.  8  by  admittance  functions  A  defined  so  that  the 
normal  components  of  the  fluctuations  are  proportional  to  the  pressure 
fluctuations , 


M'  •  d  =  ATI 


(9) 


-  +  iA^  =  lAje1^  so  that 


In  general,  A  is  a  complex  function,  A  =  A 
the  real  part  gives  that  part  of  the  normal  velocity  fluctuation  ur  which 

is  in  phase  with  the  pressure:  tu. '  =  j  p ' .  But  the  rate  at  which 

\  VP'  ,  ,  fT)  a,  ,,2 

work  is  done  on  or  by  the  waves  at  the  end  surfaces  is  ur  p  *  A'  -^s(p  )  • 

YP 


Consequently,  it  is  the  real  part  of  the  admittance  function  which  will 
be  related  to  the  damping  of  the  waves. 


For  the  cases  treated  here  (and  usually  in  rocket  chambers),  the 
boundaries  and  the  mean  flow  cause  only  small  perturbations  from  the 
classical  closed  chamber  modes  which  satisfy  the  equations 

+  h,\  - 0  (10) 


*  •  *  0  {'  ■  V  (11) 

where  R  is  the  radius  of  the  cylindrical  chamber.  For  these  modes  in  a 
closed  chamber,  the  normal  gradient  of  pressure  vanishes  at  the  ends, 

since  Vp'  =  ypVT]  =  here.  But  by  Eq.  3  with  M  =  0  and  m'  varying 

harmonically 

(bu'/dt  =  ikM7),  M'  =  -VTl/ik,  or  u'  =  -  =  -  ^7^ 

Hence  since  n  •  VYN  vanishes  on  the  boundary,  so  does  the  normal  component 
of  the  velocity  fluctuation.  In  the  case  of  purely  axial  oscillations, 

Yn  =  cos  (kji?)  with  k*  ■  At,  A  -  1,  2,  •  •  ••  The  idea  of  the  following 
computation  is  that  even  with  a  nozzle  or  a  porous  head  end,  the  mode 
shapes  are  not  much  changed  from  these  classical  vibrations;  this  is  true 
if  the  mean  flow  Mach  number  is  small  and  if  the  admittance  functions  are 
small.  It  is  not  generally  true  if  the  oscillations  are  driven  at  the 
head  end,  for  the  velocity  fluctuation  at  that  end  may  be  given  any  value 
and,  correspondingly,  the  frequency  may  be  changed  at  will. 

The  computation  here  is  directed  to  determining  the  corrections  to 
k_£  =  Att  due  to  the  mean  flow  and  damping  within  the  chamber  or  at  the 
boundaries.  It  is  a  simple  matter  to  compute  the  first-order  approximation 
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to  k  (the  mean- flow  Mach  mnber  is  the  small  parameter).  Multiply  Eq.  5 
by  Yjj,  Eq.  10  by  'll  and  subtract  the  equations;  then  Integrate  over  the 
volume  of  the  chamber: 

Mis*  -  + if  -  )JF v*  -  jflT^ 

The  volume  integral  on  the  left-hand  side  may  be  transformed  by  use  of 
Green's  theorem,  and  Eq.  6  and  11  used  for  the  surface  conditions.  One 
finds 


v» 


V* 


• SI *v) 


There  are  really  two  surface  Integrals  Involving  f— one  over  the  plane 
z  =  0  and  one  over  z  =  1. 

To  obtain  a  first-order  approximation  for  the  right-hand  side,  the 
zeroth  order  approximation  for  Tj ,  11  «  Yjj  is  used.  In  the  case  of  purely 
axial  oscillations,  Yjj  =  cos  (k/z)  with  k£  =  /Jt,  an^  if  one  assumes 

M  to  be  axial  only,  ft  =  M$z,  then  h  and  f  are  approximately  (i.e.  ,  to 
order  M) 


where  the  zeroth  order  approximation,  derived  from  Eq.  3,  has  been  used 
used  for  M  :  M7  (i/kjj)(df|/dz).  (Vote  that  the  exponential  time  factor 

exp  (iwt)  =  exp  [i(ufL/a)  •  (at/Li]  is  exp  (ikfc)  in  dimensionless  form, 
and  to  zeroth  order  k  =»  kg  =  tv.)  The  volume  integral  in  Eq.  12  may  now 
be  simplified  by  carrying  out  the  integral  over  z  with  some  use  of  the 
fact  that  d2Y//dz2  =»  -k^Y/.  One  finds 


The  surface  integral,  over  z  ■  0  and  z  ■  1,  gives 
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JJ'S*  ‘  (Jf^U  * 


lk. 


[(If  + (/"/*)  J  -  [ir»>I 


After  these  results  have  been  used  in  Eq.  22  the  real  and  imaginary  parts 
of  k  =  0  +  iA  are,  to  first  order  in  fi,* 


’i8£\R) 


y  2  dS  | 

*  V 


(15) 


*  »/  %  ♦JR) »/  ^ 


(16) 


in  vhich  Sc  is  the  cross-sectional  area  of  the  chamber  and  A  *  A  +  1A 
in  complex  form.  The  admittance  function  for  the  nozzle  end  is  An  and  at 
the  head  end  it  is  A0. 

It  is  now  apparent  why  the  calculations  were  not  immediately  reduced 
to  one-dimensional  form.  If  that  step  had  been  taken,  then  the  important 
(and  obvious)  result  that  the  surface  integrals  must  be  weighted  by  the 
appropriate  areas  would  not  have  appeared.  Thus,  for  example,  if  the 
nozzle  entrance  area  is  S_,  the  damping  term  in  Eq.  16  is  An'f’^S_/sc, 
since  =  1  at  z  =  1. 


The  formulas  15  and  1 6  are  valid  for  the  axial  modes  with  one- 
dimensional  mean  flow,  subject  to  the  restriction  that  the  Mach  number  of 
the  mean  flow  must  be  small.  Hence,  care  must  be  taken  when  these  are 
used  for  a  rocket  with  an  abrupt  change  of  cross  section,  from  the  port 
area  to  a  sonic  nozzle  as  sketched  in  Fig.  1.  There  are  evidently  regions 
of  rapid  variations  in  mean-flow  properties  near  the  porous  head  plate  and 
near  the  nozzle.  Average  streamlines  have  been  sketched.  The  expressions 
15  and  1 6  are  valid  only  in  the  region  between  the  planes  at  a  and  b, 
between  which  the  Mach  number  is  approximately  constant,  equal  to  Mg  when 
there  is  no  burning. 


It  appears  that  the  region  near  the; porous  plate  may  be  so  thin  that 
as  an  approximation,  position  a  is  in  fact  at  the  origin,  and  the 


* 


Note  that 


Y/dz  =  £  for  /  /  0. 
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admittance  function  is  then  non-zero  over  the  porous  area  (Sp)  at  most 
(the  velocity  fluctuation  is  elsewhere  zero  on  z  »  0).  But  yet  the 
Influence  of  the  mean  flow  may  be  represents*  by  the  average  chamber  Mach 
number  acting  over  the  entire  chamber  area.  Mote  that  the  volvmm  Integral 
of  h  and  the  part  of  f  depending  on  Mach  number  cancel,  which  happens  only 
if  M  Is  constant.  This  is  merely  a  statement  of  the  fact  that  if  the  flow 
is  uniform,  the  rate  of  convection  of  wave  energy  Into  the  region  bounded 
by  two  planes  normal  to  the  flow  is  obviously  equal  to  the  rate  of  con¬ 
vection  out.  Then,  If  the  admittance  functions  are  constant  over  the 
corresponding  areas,  Eq.  15  and  l6  become 


To  first  order,  then 


k~kA 


The  effects  at  the  exhaust  end.  Including  rapid  acceleration  to  the 
sonic  throat  and  eddies  of  some  sort  which  must  exist  in  the  corners  of 
the  chamber,  have  been  Included  In  Aq.  It  is  obviously  not  reasonable  to 
calculate  this  quantity  which  Is  more  than  Just  the  nozzle  admittance 
function.  But,  on  the  other  hand,  if  one  tries  to  extend-  the  region  In 
which  the  standing  waves  exist  (i.e.,  b  ""  l),  then  tjhe  equations  for  the 

waves  must  be  seriously  complicated,  both  because  |M|  is  not  small  and 
because  M  is  no  longer  in  the  axial  direction.  The  analysis  by  Cantrell, 
McClure ,  and  Hart  (Ref.  6)  is  of  course  subject  to  the  same  restrictions. 
If  the  nozzle  convergent  fairs  more  or  less  smoothly  to  the  chamber,  then 
Ajj  is  indeed  almost  entirely  the  nozzle  admittance  function,  which  can  be 
calculated  independently. 


It  is  these  difficulties  which  lead  one  to  vise  the  analysis  as  a 
guide  to  interpretation  of  experimental  results,  rather  than  as  a  means 
of  thorough  computation.  According  to  Eq.  l6  and  17,  the  amplitude  of  a 
wave  varies  with  time  as 


In  terms  of  dimensional  quantities,  the  attenuation  constant  a  is 
therefore 


t&C&tS&K 


HWC  IP  451*4 


“  - 1  (AJ 


(r)  i 


and  the  angular  frequency  of  flteely  decaying  waves  Is 
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Equations  19  and  20  are  applicable  to,  for  example,  experiments  Involving 
the  explosion  of  a  small  charge  and  measurement  of  the  decay  rate  and 
frequency  of  subsequent  waves.  They  are  valid,  in  fact,  for  a  port  of  any 
shape  so  long  as  It  Is  straight.  If  the  Mach  number  Is  not  uniform,  some 
of  the  integrals  In  Eq.  15  and  16  contribute  additional  terms. 


IV.  WAVES  WITH  UNIFOFM  MEAN  FLOW 


has 


For  studying  the  frequency  response  when  waves  are  driven,  it  is  best 
to  incorporate  at  the  very  beginning  the  fact  that  the  mean  flow  is  uni¬ 
form  (although  as  in  the  previous  section,  the  cross  section  of  the 
chamber  need  not  be  circular) .  The  governing  first-order  equations  for 
M  ,  the  axial  Mach  number  fluctuation,  and  T|,  are 

<a> 

<-> 


3 


in  which  M.  is  the  Mach  number  of  flow  in  the  chamber.  The  case  of 
transient  decay  of  waves,  i.e.,  the  transient  response  of  the  chamber, 
has  already  been  treated,  with  the  results  given  in  Eq.  19  and  20.  In 
this  section,  the  analysis  of  the  frequency  response  is  begun. 

By  analogy  with  the  familiar  ideas  of  "lumped"  electrical  circuits 
and  mechanical  systems,  one  expects  the  width  of  the  peak  in  the  fre¬ 
quency  response  to  be  proportional  to  the  attenuation  constant,  Eq.  19, 
associated  with  the  transient  response.  Some  experiments  (Ref.  5  for 
example )  have  in  fact  shown  quite  good  agreement,  but  it  is  not  obvious 
why  this  should  be  so.  Indeed,  the  calculations  here  indicate  that  the 
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agreement  may  imply  that  the  admittance  function  for  the  porous  plate 
(and  hence  the  associated  attenuation)  Is  negligible  compared  to  the 
effects  of  the  exhaust  end  of  the  chamber.  This  was  Intended  to  be  the 
case  in  the  experiments  (Ref.  4). 

When  standing  waves  are  sustained  by  external  means,  the  analysis  is 
often  instructively  cast  In  terms  of  traveling  waves  (as,  for  exaaple,  in 
the  problem  of  a  transmission  line  which  the  present  problem  closely 
resembles ) .  If  one  performs  the  coordinate  transformation 

T  =  t 

C  =  z  -  M  t 
c 

then  Eq.  21  and  22  become 


+  an  n 

dT)  2m'  n 

3r  +3T  0 


which  imply  the  wave  equations 


(fr '  fr)  t ) 


Obviously  then,  one  has  two  solutions,  corresponding  to  waves  traveling 
to  the  left  and  to  the  right.  Alternatively,  the  wave  equations  in  t, 
z  coordinates  are 


-os 


#  #  1/m 


cE3t  '  ^2 


](;') 


In  any  case,  the  complete  solutions  can  be  written 


11  -  eJ 


„/  lPt 
M  “  e 


[v'1  rV  ♦ -1  rV] 

[„  -i  n  .  ..  i  n  1 

M+  e  t-rr z  +  M- e  rnr  ZJ 


Rote  that  the  frequencies  of  the  two  waves  are  not  the  same: 

“0/(1  +  M  )  for  the  wave  to  the  right,  and  -  fi/(l  -  M  )  for  the 
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wave  to  the  left.  This  is  the  "Doppler  effect"  often  mentioned  in  the 
literature  of  acoustics  applied  to  rockets,  but  not  usually  shown.  A 
wave  of  frequency  0  in  the  laboratory  and  traveling  in  the  direction  of 
the  moving  stream  has  a  lower  frequency  with  respect  to  the  stream. 
Substitution  of  these  expressions  into  one  of  the  first-order  Eq.  21  or 
22  shows  that 

M+  =  \ 

M_  =  -T|_  (27) 

Now  in  the  case  of  the  cold-flow  rocket  experiments,  the  waves  sire 
driven  at  the  head  end.  A  solution  corresponding  to  the  standing  wave 
existing  at  any  frequency  is  obtained  by  superposing  the  correct  amount 
of  +  and  -  waves  to  satisfy  the  boundary  condition  at  the  exhaust  end, 
assumed  here  to  be  at  a  plane  as  suggested  by  Fig.  1  and  associated 
remarks.  One  can  imagine  that  a  wave  traveling  to  the  right  is  partly 
reflected  and  partly  transmitted;  the  amplitude  and  phase  of  the  reflected 
wave  is  so  determined  as  to  have  a  standing  wave  observed  in  the  labora¬ 
tory.  The  ratio  of  amplitudes  of  the  impinging  and  reflected  waves  is 
usually  called  the  "standing  wave  ratio . "  In  accord  with  previous 
remarks,  let  all  effects  at  the  exhaust  end  be  represented  by  the 

admittance  function  A  ,  so  that  at  z  =  1 

n’ 

M1/1!  =  A  (28) 

n 

The  plane  z  =  1  is  actually  slightly  displaced  from  the  end  plate--it  is 
labeled  b  in  Fig.  1.  Note  again  that  includes  more  than  simply  the 
nozzle  admittance  function.  Hence  Eq.  25-27  imply 

\  e-1"  -  \  elS 
\  e'^  *  1)_  e1" 

where  0  =  Q/(l  -  Mc^)  and  use  has  been  made  of  the  fact  that 

t— =  m  n  +  n,  ■  v-  =  -m  5  +  q 

l-M  c  '  1  +  M  c 

c  c 

Hence  the  complex  amplitude  of  the  wave  traveling  to  the  left  is 


1  -  A 

•n  „  n  -2i0 

\  - — e 
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T!  -  T| 
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which  can  also  be  written  as 

an 


♦  m<Hc*  -  1)  S{,  -  1)  -  .In  ft(,  -  l)] 


M 


pz-r  .int  +  -  l>  [a,  co.  n(i  .  1)  -  1  sin  n(i  - 1)] 


Now  define  4  by  tanh  4  *  and 


2TK 


(A  +  1}  cosh  4 


Si°t  +  ^(“c*  "  1)  cosh  [in(z  -  1)  +  4]  (29) 


M' 


Ta 


— -  iTTS.h  if' eint  *  in(McZ  '  l)  ■*»»  (30) 


These  results  will  be  used  to  describe  Interpretation  of  measurements  of 
the  steady«8tate  frequency  response  and  of  the  amplitude  of  oscillations 
along  the  chamber  axis. 


V.  FREQUENCY  RESPONSE 

The  measurements  have  been  made  of  pressure  at  the  head  end.  It  may 
be  supposed  that  the  oscillations  are  driven  by  fluctuations  In  velocity, 
or  Mach  number,  at  the  head  end.  Thus,  let  M  *  Mq'  e^  at  z  ■  0,  and 
Xq.  29  gives 


2\ 


cosh  ^(A^  +  1) 


-M  '  c1" 

_ o _ 

■ 

slnh  (4  -  ifl) 
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Then  the  measured  frequency  response  is  represented  by 
T|  cosh  [ifl(z  -  l)  + 

-  ifi) 

In  particular,  at  z  =  0, 


(31) 


cosh 

sinh  -  in) 


which  can  be  expanded  to 


1  +  iA  tan  n 
n 

A  +  i  tan  ft 


(32) 


It  is  particularly  important  to  note  that  since  the  waves  are  driven 
by  perturbations  at  the  porous  plate,  the  influence  of  that  region  (con¬ 
tained  in  Plq  defined  earlier)  cannot  appear:  they  are  combined  with  the 
external  driving  mechanism  to  provide,  for  a  given  frequency  and  magni¬ 
tude  of  pressure  fluctuation,  just  the  value  of  Mq '  required  to  maintain 
a  steady  oscillation.  Hence  measurements  made  with  a  transducer  in  the 
porous  plate  can  yield  no  information  about  the  admittance  function  to 
be  associated  with  the  porous  plate  when  waves  are  decaying  freely  within 
the  chamber.  Let  =  u  +  i0: 


\  O  /3=0 


(1-9  tan  fi, 

'-1 — rm* 


.  1  x  “ 
+  tan  n) 


and 


|Z|  = 


Im. 


fW 

’  D 


where 


N  =  X  tan^  Q  -  20  tan  fi  +  1 
D  =  tan^  n  +  20  tan  0  +  X 
X  =  u2  +  02 


(33) 


(3«0 

(35) 

(36) 
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The  function  |z(Q)|/Jz(Q)  jn-p  la  plotted  in  Fig.  2 -6  for  several 
representative  values  of  |i  and  6 •  Peaks  occur  at  frequencies  fi=  Cl 
which,  when  both  4  and  3  are  small,  are  close  to  the  frequencies  of  the 
natural  model  for  a  closed  chamber  (Op  **  in)  and  close  to  the  frequencies 
of  transient  oscillations  (Eq.  20)  but  not  necessarily  equal  to  either. 
The  peaks  are  determined  from  the  condition 


d  Z 


=  0 


(37) 


To  simplify  calculations,  and  also  to  avoid  introducing  further  unknown 
quantities,  suppose  that  both  9  and  4  are  slowly  varying  functions  of  0 
near  the  peak,  then  both  d.9/dP  and  du/dO  will  be  set  equal  to  zero.  One 
finds  then  that  C2p  is  the  solution  to 


tan^  0 


Thus 


1  -  X 

tan  %  55  *20“ 


(38) 


+  4 


(39) 


The  sign  on  the  radical  cannot  be  fixed  uniquely.  In  the  limit  of 
small  8  and  4  the  negative  sign  Is  the  correct  choice,  for  then  Eq.  39 
gives 


8 


tan  ft  -  .  v 
p  1  -  X 


-8 


and  hence 


Q  «*  In  -  8 
P 


(40) 


which  Is  exactly  the  correction,  due  to  the  nozzle,  shown  in  Eq.  20  for 
the  frequency  of  the  transient  response.  On  the  other  hand,  peaks  do 
exist  for  0,  far  removed  from  j8tt-  a  particularly  interesting  limit  is 
found  by  taking  the  positive  sign  in  Eq.  39,  with  9  small,  for  which 
Eq.  39  becomes 


tan 


1  -  X  1  -  4 


(41) 


Let  1  -  \F  be  non-zero,  and  for  9  -*  0  the  peak  frequencies  fall  in  the 
neighborhood  of  multiples  of  tt/2,  as  shown  later  by  the  numerical 
results . 
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FIG.  5.  Dimensionless  Impedance  at  Head  End  as  a  Function  of 
Frequency  0  =  0.4. 
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Of  greater  practical  interest  is  the  width  of  the  peak,  for  one 
expects  by  analogy  with  electrical  circuits  and  simple  mechanical  systems 
that  this  should  be  related  to  the  damping  of  the  system.  Experiments 
(Ref.  3  and  4)  have  shown  that  the  peak  is  oft®n  narrow;  the  width  at  the 
"half-power  points,"  i.e.,  the  frequencies  at  which  | Ti/Mq '  \  is  down  by 
l//2  will  be  denoted  by  2e.  For  0  and  n  small,  the  frequencies  at  the 
half -power  points  are  therefore  approximately  in  +  S  _+  e.  For  this 
limiting  case,  the  peak  value  of  I Z 1  is 


and  one  finds 


(n,e  «  1) 


Hence,  under  these  conditions,  the  width  of  the  peak  (2e  2 u)  is 

twice  the  attenuation  constant  for  the  transient  oscillations  (Eq.  19) 
providing  the  porous  plate  does  not  cause  much  damping  of  the  transient 
motions  (i.e.,  AgC1")  =  0).  It  follows  that  in  this  case,  the  damping  at 
the  exhaust  end  may  indeed  be  found  from  Q  =  ftp/2e;  this  is  valid,  then, 
when  the  response  has  a  narrow  peak,  associated  with  relatively  light 
damping.  The  attenuation  constant  is 

n 

a=_E  =  £_E  =  S  (2€)  (k2) 

2 Q  L  2Q  L  K  J  K 


It  appears  therefore  that  if  agreement  is  found  between  the  value  given 
by  Eq.  4l  (i.e.,  by  measurement  of  the  frequency  response)  and  the  value 
given  by  Eq.  l8  (i.e.,  by  measurement  of  the  transient  response)  then  the 
real  part  of  Aq  must  be  negligibly  small.  That  is,  assuming  always  a 
lightly  damped  system,  so  that  the  approximations  leading  to  Eq.  40  and 
4l  are  valid,  one  would  conclude  that  the  porous  plate  contributes 
insignificant  damping. 
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Not  all  experiments  yield  narrow  peaks  in  the  frequency  response  and 
in  those  cases  the  approximate  calculation  used  above  is  not  valid.  One 
cannot  compute  the  transient  attenuation  constant  from  the  simple  for¬ 
mula  42  involving  Q.  It  is  necessary  to  use  some  sort  of  graphical  pro¬ 
cedure  such  as  the  following.  The  values  of  |z|  at  the  half -power  points 
are  still  given  by 

|z4  "k  |Z-L*  *72  Vi£ 

where  Np,  Dp  are  the  values  of  N,  D  for  0  =  fL.  This  equation,  with  N  and 
D  given  by  Eq.  34  and  35,  leads  to  a  quadratic  equation  for  tan  0,the 
tangent  of  the  normalized  frequency  at  the  half-power  points: 


tan  0  +  tan 


'-20 (D 

n  — -k 

2XD 

l 


N 


[D  -  XN  /2*1 

-2 _ EL 

XD  -  N  /2 

.  P  P  J 


=  0 


(43) 


From  Eq.  43  one  can  find  tan  0  and  hence  fi  at  the  upper  and  lower  half¬ 
power  points  as  functions  of  p  and  0.  The  results  are  conveniently- 
displayed  in  the  forms  shown  in  Fig.  7“ 10. 

Those  graphs  can  be  used  to  determine  both  the  real  and  imaginary 
parts  of  the  admittance  function  at  the  exhaust  end  by  measurement  of  the 
frequency  of  the  peak  and  the  width  of  the  response  function.  The  same 
information  has  been  cross-plotted  in  several  ways;  it  is  not  clear  at 
this  time  just  which  form  may  be  most  convenient  for  use  in  the 
interpretation  of  data. 

Note  that  in  Fig.  2-6,  the  impedance  has  been  normalized  with  respect 
to  its  value  at  zero  frequency.  Those  charts  also  6how  a  feature  remarked 
upon  earlier  in  connection  with  Eq.  39,  that  the  frequencies  of  the  peaks 
need  not  be  close  to  the  values  for  a  closed  chamber  in  the  absence  of 
flow.  However,  it  appears  that  in  practice,  both  p  and  0  are  small  so 
that  in  fact  Op  lies  near  In,  as  shown  by  Eq.  40. 

In  Fig.  7,  there  are  several  limiting  values  which  should  be  men¬ 
tioned.  The  case  9=0  yields  a  resonant  frequency  of  in  if  p  <  1  and 
irr/2  if  4  >  1.  This  behavior  may  also  be  seen  in  Fig.  2,  and  follows 
from  Eq.  40  and  4l,  which  show  that  the  switch  between  values  occurs  at 
p  =  1.  The  limiting  case  9  -*  00  can  also  be  extracted  from  Eq.  4l. 

A  quantity  called  "maximum  half-power  bandwidth, "  equal  to  tt  appears 
in  Fig.  8-10.  This  is  also  a  limiting  case,  which  arises  because  the 
impedance  Z  is  precisely  periodic  in  with  period  tt.  That  the  period 
i6  tt  can  be  shown  easily  by  substituting  0  +  tt  for  fi  in  Eq.  33*  This  is 
exactly  true  only  if  p  and  9  are  both  taken  to  be  independent  of  frequency. 
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FIG.  8.  Half -Power  Bandwidth  of  the  Impedance  at  the  Head  End. 
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Half -Power  Bandwidth,  &2  -  &  l 


FIG.  9.  Frequency  at  the  Maximum  Value  of  the  Impedance  at  the  Head 
End  as  a  Function  of  Half-Power  Bandwidth.  ( —  Constant  u, 

—  Constant  8.) 
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an  assumption  used  In  all  the  graphical  results.  However,  this  means 
that  the  distance  between  peaks  Is  tr  always,  and  therefore  the  marl  wan 
bandwidth  la  also  ft,  being  tt/2  on  each  side  of  a  given  peak. 


VI.  AMPLITUIE  OF  06CILLATXGR8  ALCRO  THE  CHAMBER; 
8TAHBH0  HAVE  RATIO 


If  provision  is  made  for  measuring  pressure  fluctuations  as  a  function 
of  position  along  the  axis,  the  real  and  Imaginary  parts  of  Aq  may  be 
determined  much  more  readily.  Equation  31  gives 


An  + 


| | cosh  Cifi(.  -  1)  +  4]| 

jeosh  [4  +  ifi(z  -  l) 3 J 
Jtosh^a  -  sin2  [p  -  5(z  -  l)] 


where 


4  -  a  +  ip  -  tanh-1(-An)  ■  tanh^C-u  -  18) 


This  last  equation  can  be  split  Into  real  and  Imaginary  parts: 


..ilniki4^4] 

4  La  - «)  *  rj 


e-5t“'1 1 .  ih? 


Equation  44  shows  that  |  T)  |  has  a  maximum  value  (for  a  given  frequency  and 
T)+,  or  Mo#)  when  0  -  -  l)  =  0,  tt,  2tt,  .  .  .,  that  is,  when 

6(1  -  z)  =  nn  -  0 


or,  in  terms  of  the  real  frequency  f  *  ^  (l  -  Mc2)  ^  cps,  and  real 
distance  x  *  zL, 


a(l  -  M^) 


(L  -  x)  *»  nn  -  0 
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Thus,  the  first  maximum  occurs  at  a  distance 


(L  -  x)  .  (-|S)(1  -  M02)  ±  (48) 

from  the  exhaust  if  0  <  0,  and  at  a  distance 

(L  -  x)  =  (tt  -  p)(l  -  Me2)  ^  (49) 

if  0  <  0  <  tt  (X  ■  a/f ) . 

Moreover,  |T)j  has  a  minimum  value  when 

- —  -  (L  -  x)  =.  (n  +  i)  n  .  0  (50) 

a(l  -  M.2)  2 

i.e.,  at  a  distance 

(L  -  x)  .  -  p)(l  -  Mc2)  ^  (51) 

from  the  exhaust  end,  if  0  <  tt/2,  or 

(L  -  x)  -  (2J  -  P)(l  -  Mc2)  ^  (52) 


if  0  <  3tt/2.  The  ratio  of  the  maximum  to  minimum  values  of  T)  is 


UL* 


cosh  a 

•/coeh^  a  -  1 


tanh  a 


(53) 


Evidently,  then,  measurement  of  the  position  of  a  maximum  or  minimum 
of  |T|  |  and  the  ratio  |T)|max/l^lmin  18  sufficient  to  determine  a  and  0 
from  Eq.  48  and  49,  or  51  and  52.  Finally,  u  and  9  can  be  computed  from 
Eq.  46  and  47.  An  alternative  approach  is  to  use  a  graphical  represen¬ 
tation  such  as  the  Smith  chart  commonly  employed  in  calculations  for 
transmission  lines.  (See  also  Chapter  VI  of  Ref.  7>)  If  these  results 
are  combined  with  measurement  of  the  attenuation  constant  of  freely 
decaying  waves,  then  the  real  parts  of  the  admittance  functions  for  both 
ends  of  the  chamber  will  be  known. 


NWC  TP  4544 


VII.  CONCLUDING  REMARKS 


It  must  be  emphasized  that  the  calculations  covered  in  this  report 
are  intended  to  be  coupled  closely  with  experimental  work  in  an  effort 
to  determine  the  acoustic  losses  associated  vlth  a  rocket  chamber. 

Although  some  contributions  to  the  dissipation  of  energy  can  be  computed 
quite  accurately,  such  as  that  due  to  viscous  shear  at  the  walls,  others 
which  have  been  noted  earlier,  are  simply  beyond  analysis.  Moreover,  the 
techniques  discussed  here  can  probably  be  extended  to  awkward  geometries. 

The  principal  results  of  the  present  analysis  are  the  following: 

1.  The  computations  in  section  III  give  simple  formulas  for  the 
frequency  and  attenuation  constant  for  each  natural  mode  (Eq.  17  and  18) 
when  the  influences  of  the  ends  can  be  lumped  into  simple  admittance 
functions.  This  implies  that  the  regions  of  transition  between  the  uni¬ 
form  flow  in  the  chamber  and  the  end  surfaces  must  be  short  compared  to 
the  chamber  length;  or,  better,  that  Ajj  and  Aq  must  be  small. 

2.  The  frequencies  of  the  natural  modes  are  not  in  general  equal  to 
the  frequencies  at  which  the  steady-state  frequency  response  has  peaks. 

The  former  occurs  at  normalized  values  in  -  (A^\rJ  +  n)  according  to 
Eq.  17.  The  latter  occur  when  Eq.  53  has  maxima.  For  a  given  mode,  the 
two  values  are  approximately  equal  when  is  small  and  Aq  is  negligible. 

3.  The  width  of  the  peak  in  the  frequency  response  is  sinply  related 
to  the  transient  attenuation  constant  only  when  A^  is  small  and  Aq  is 
negligible.  If  A^  is  not  small  (i.e.,  the  peak  is  broad,  low  Q),  then 
the  width  of  the  peak  is  related  to  both  the  real  and  imaginary  parts  of 
An,  but  interpretation  is  best  done  graphically  by  making  use  of  Fig.  2-10. 

4.  Measurement  of  the  positions  of  maxima  and  minima,  and  the  ratio 
of  maximum  to  minimum  ("standing  wave  ratio")  will  give  the  real  and 
imaginary  parts  of  Ajj  directly.  If  the  measurements  can  be  made  without 
great  difficulty,  then  this  seems  to  be  a  better  means  of  determining  the 
losses  associated  with  the  exhaust  end  of  the  chamber. 


5.  The  experimental  results  so  far  reported  indicate  that  the  losses, 
and  the  shift  of  frequency,  due  to  the  exhaust  flow  are  very  much  larger 
than  those  due  to  the  porous  plate.  That  is,  under  most  conditions,  the 
porous  plate  seems  to  offer  "infinite  impedance."  However,  the  validity 
of  these  results  is  questionable  at  the  present  time  owing  to  some 
internal  inconsistencies;  thus  the  data  have  not  been  examined  here. 

Determination  of  the  acoustic  losses  in  a  rocket  by  making  measure¬ 
ments  in  a  laboratory  scale  model,  without  burning,  appears  at  the  present 
time  to  be  a  very  promising  aid  to  design.  However,  the  experimental  work 
has  not  yet  been  sufficiently  developed  to  make  extensive  use  of  the 
calculations  discussed  here. 
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